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ABSTRACT 

We present the recent discovery of a new subdwarf B variable (sdBV), with an exceptionally low surface 
gravity. Our spectroscopy places it at r eff = 32100+1000, logg = 5. 15+0.20, and log (NhJNh) = -2.8+0.2. With 
a magnitude of B = 12.0, it is the second brightest V361 Hya star ever found. Photometry from three different 
observatories reveals a temporal spectrum with eleven clearly detected periods in the range 376 to 566 s, and 
at least five more close to our detection limit. These periods are unusually long for the V361 Hya class of 
short-period sdBV pulsators, but not unreasonable for p- and g-modes close to the radial fundamental, given 
its low surface gravity. Of the ~50 short-period sdB pulsators known to date, only a single one has been found 
to have comparable spectroscopic parameters to GALEX J201337.6+092801(J20136+0928, for short). This 
is the enigmatic high-amplitude pulsator V338 Ser, and we conclude that J20136+0928 is the second example 
of this rare subclass of sdB pulsators located well above the canonical extreme horizontal branch in the HR 
diagram. 

Subject headings: subdwarfs - stars: early-type - stars: variables: general - stars: individual: GALEX 
J201337.6+092801 
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1. INTRODUCTION 

The subdwarf B (sdB) stars are generally known to be 
core-helium-burning stars with inert hydrogen-dominated en- 
velopes of very low mass. This composition places them on 
the hot extension of the horizontal bran ch, the extreme hori- 
zontal branch (EHB: lHeberi[l9 86, 2009). Their low-mass en- 
velopes prevent them from reaching a second giant stage, and 
instead, after their core helium is exhausted, they expand only 
briefly, never exceeding more than about a third of a solar 
radius, before contracting and passing on towards higher ef- 
fective temperatures and gravities. Thus, an sdB star evolves 
into the hotter sdO population before reaching d egeneracy and 
the as sociated white dwarf (WD) cooling track (iDorman et al.l 
fl993h . 

While the future evolution of sdB stars is unproblematic, 
exactly how they lose their envelopes during their first red 
giant (RGB) stage is more complicated. If the star has a close 
low-mass companion, mass transfer on the tip of the RGB 
will quickly pull it into the envelope. As friction transfers 
angular momentum from the orbit to the envelope, it will spin 
up and be ejected from the system. If the core has become 
massive enough for helium to ignite (~O.47M ), the result 
will be a very close sdB+dM system. These binaries have 
spectacular reflection effects, and a total of 14 such object s are 
described in the literatu re; 13 are listed inlFor et al.l(l2010l) and 
the most recent one in IQstensen et al.l d2010al) . While such 
sdB+dM binaries make out only a tiny fraction of the total 
sdB population, cases where the companion responsible for 
ejecting the en velope is a low-mass WD, appear to be much 
more common. iMaxted et al.l (12001b found that 21 of 36 stars 
in their radial velocity study are in short period systems, most 



likely with WD companions. 

If the companion is more massive than the subdwarf (at 
least at the end of mass transfer), the orbit will have ex- 
panded to more than a hundred days. Such orbits are hard 
to measure, but the companion is easily detectable spectro- 
scopically or from infrared excess. iNapiwotzki et al.l (120041) 
found that more than a third of their sdB sample show the 
spectroscopic signature of main sequence companions, while 
iReed & Stienind ( 120041) . using 2MASS photometry, inferred 
that about half of the sdBs in the field are likely to be of this 
type. Thus, the vast majority of sdB stars in the field must be 
in binary systems. However, a few pulsating sdBs have been 
studied in such detail that any companion close enough to 
have interacted during the RGB stage can be ruled out. A pos- 
sible explanation for such single sdB stars can be the merger 
of two low-mass WDs. These three scenarios are all inves- 
tigated in de tail in the binary population study of lHan et al.l 
d200Hl2003l) . They find that while the first two scenarios must 
have sdBs with masses close to the He-flash mass of 0.47 M , 
the merger scenario can produce sdBs with a broad distribu- 
tion of masses ranging from 0.4 to 0.7 M . 

While most sdB stars form a band in the r e ff-log g diagram 
that is consistent with a mass close to the canonical 0.47 M , 
and various envelope thicknesses up to the H-burning limit, a 
few exceptions are found. One such is V338 Ser, discove red 
to be a high amplitude sdB pulsator bv lKoen et alj (119981). A 
multi -site campaign devoted to this pulsator dKilkennv et al.l 
[19991) found more than fifty pulsation frequencies in the 
range 365-529 s, which now defines the long-period edge for 
the V361Hya group of short-period sdB variables (sdBVs). 
Time-resolved spe ctroscopy of V33 8Ser was analysed by 
iTillichet all (120071) . and most recentlv lOstensen et all ( 1201 Pel ) 
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FIG. 1. — Light curves of J20136+0928. The top panel represents the discovery light curve. In the online version the different data sets are shown in different 
colors (Mercator; red, Baker; green, SPM; blue). 



performed high-resolution time-resolved spectros copy. Its 
unusu ally low gravity was discussed briefly by I0stensen1 
proposing that it may be a core-He-burning EHB star 
with a mass larger than the canonical, rat her than a shell-He - 
burning post-EHB star as sugg ested by iKoen et al.l (119981) . 
Asteroseismic solu tions supporting this hy pothesis were re- 
cently presented bv lvan Grootel et all (120101) . implying a mass 
of 0.76 M , and an envelope mass fraction of 0.2%. 

Since V338 Ser remained a unique object in spite of a 
decade of searches for V361 Hya pulsators, we have recently 
been paying particular attention to candidates in this low- 
gravity region. Here we present the discovery of pulsations 
in J20 136+0928 (Figure [TJ, which has an even lower surface 
gravity than V338 Ser, and the lowest of any sdB pulsator 
found to date. 

2. TARGET SELECTION 

Since 2008 we have observed UV-excess targets selected 
based on UV photo metry from the Gala xy Evolution Explorer 
(GALEX) satellite (iMartin et al.l 120051) . An estimate based 
on targets with NUV < 14.0 and FUV-B < +0.5 is that only 
about two-thirds were already cataloged. UV-excess objects 
in this cut include a host of interesting compact stars. All 
types of pulsating WDs and hot subdwarf stars are included, 
as are most types of cataclysmic variables and planetary neb- 
ula nuclei. Our first cut based on GALEX data release 4 (GR4), 
contained 649 objects of which only about 400 had even ap- 
proximate classes in the literature. We started observing tar- 
gets from this list as a backup programme for various observ- 
ing runs, and have by now mostly completed the sample. Full 
details of this survey will be give n in a forthcoming p aper. 
A similar survey is described by IVennes et al.l (1201 ll) . As 
we have accumulated spectroscopic observations, new sdBV 
candidates have entered our observing lists for photometric 
follow-up, and new discoveries are emerging. The first result 



was the di scovery of pulsatio ns in J08069+1527, recently pre- 
sented bv lBaran et al.l d201 ll) . and establishing it as a being a 
high-amplitude member of the rare DW Lyn (hybrid) type of 
sdBVs. 

The object presented here appears in GALEX data release 6 
(GR6) as GALEX J201337.6+092801, with UV magnitudes 
FUV= 12.263(3), NUV= 13.079(3). Photographic magni- 
tudes for this star from the NOMAD database are B = 12.01, 
V= 12.22, tf= 12.34, and from the 2MASS near-IR survey; 
/= 13.00, H= 13.11, K= 13.23. The magnitudes are all per- 
fectly consistent with an sdB star that has no main sequence 
companion. J20 136+0928 also appears in the Tycho-2 catalog 
as TYC 1077-218-1, where it has B= 12.16, V= 11.78. The 
NOMAD database also lists the target with a small but signifi- 
cant proper motion (/i [V ,yUtf = 12.1,2.8 mas/yr). 

The target was observed spectroscopically during a run at 
the Isaac Newton Telescope (INT) at Observatorio Roque de 
los Muchachos on the island of La Palma (Spain). After the 
observing run we processed all the data and extracted the 
spectra with IRAF in the regular way. We fitted the hot sub- 
dwarfs using the pure H+He LTE model grids of [Heber et al.l 
(2000), and we show the resulting fit for J20 136+0928 in Fig- 
ure |2] The formal fitting errors (shown on the figure) do not 
account for systematic effects inherent in the models, so we 
generously increase the errors by a factor five when stating 
r eff = 32 100+1000, logg = 5. 15+0.20, and log (iV H e /#£)=- 
2.8+0.2. Note that the typical shift of when going from LTE 
to NLTE models is ~+1000K, and can bec ome larger w hen 
metalicity effects are taken into account (see lHeberll2009l and 
references therein). However, the low surface gravity of this 
object is clear. 

3. PHOTOMETRIC OBSERVATIONS 

The discovery run was made on October 6, 2010 with the 
1.2-m Mercator Telescope, sited just next to the INT. Al- 
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1 J20 136 + 0928 T elf = 32095±soi 

log g = 5.150±0.040 
. I 10JE log He/H = -3.763±0.065 



FIG. 2.— LTE model fit of J20136+0928. 

though the conditions were poor, the main period of ~520 s 
was clearly seen, and recognised to be longer than normal for 
a V361Hya star. Since it was already clear from the spec- 
troscopy that the gravity of the target was unusually low, the 
connection with V338 Ser was immediately made, and we ur- 
gently proceeded to confirm the discovery. As the target was 
already setting early in the night, we pooled observing time at 
several telescopes in order to collect enough photometry for a 
frequency analysis. 

Photometric data were collected at three different sites 
within one month of the discovery run. At the Mercator 
telescope we used the Merope CCD camera, which was re- 
cently upgraded with a new E2V frame trans fer CCD with 
2048x3074 illuminated pixels (l0stensenl2OlOh . The observa- 
tions were done with a Geneva-B filter (A = 420 1+286), and 
a cadence period of ~15 sec. With the 0.4-m telescope at 
Baker Observatory (Missouri, USA) we used a Photometries 
RS-1340 CCD with a BG40 filter (,1 = 4750+1500), and a ca- 
dence of ~30s. We also obtained three consecutive nights 
with the 0.84-m telescope at San Pedro Martir (SPM) obser- 
vatory in Mexico. The CCD camera is equipped with an E2V 
4240 CCD, and we used a Bessel-fi filter (A = 4350+980 ), and 
a cadence of 38 s. 

The light curves are all shown in Figure Q] with the discov- 
ery run on top. Poor weather at Mercator left us unable to 
collect a significant amount of data at our prime site. Differ- 
ential photometry was made with respect to the sum of three 
or four of the brightest reference stars found within the field 
of view of the CCDs. Most of this reference signal comes 
from TYC 1077-402-1, about 1.7' SW of the target, which is 
comparable in brightness to J20 136+0928. The five bottom 
panels are consecutive nights, but the coverage is rather poor 
due to the target only being visible for a few hours before set- 
ting. On both October 27 and 28, there was a few hours of 
overlap between data from Baker and SPM. This allowed us 
to confirm that there is no significant difference in amplitude 
between observations made with the two telescopes, in spite 
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FIG. 3. — FT of the five consecutive nights taken together (top panel), and 
the prewhitening sequence after fitting and subtracting five, ten and sixteen 
frequencies from the original light curve. The window function is shown as 
an inset to the top panel. The arrows indicate the positions of the frequencies 
identified and removed in each step. The scale in the bottom panel has been 
expanded to show the residual level, and the line at which we stop is at 3.7 
times the mean residual level. 

of a significantly broader bandpass at Baker. We only use the 
last five nights for the frequency analysis described in the next 
section. 

4. FREQUENCY ANALYSIS 

The Fourier transform (FT) of all runs longer than 1 h is 
far too short to resolve the complexity of the frequency spec- 
trum displayed by J20 136+0928. Joining the five consecutive 
nights and taking the FT of the combined dataset produces 
the complex frequency spectrum with strong 1-day aliasing 
shown in the top panel of Figure|3] The tools used for the fre- 
quency determination were d eveloped particularly for multi- 
site data, and are described in Vuckovic et al.l (1201 ll) . We pro- 
ceed to prewhiten the most significant well-separated peaks, 
five or six at a time, until we reach a level of 0.71 mma 
(straight line in the bottom panel of the same figure). This 
level is 4.2 times the noise level when computed in the high- 
frequency region where no pulsational power is seen (al- 
though we do detect the drive frequency of some of our tele- 
scopes), and 3.7 times the mean residual level of the 1000 
to 3000 /Aiz region shown in the last panel. In each step the 
most significant peaks are used as input for a non-linear least- 
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TABLE 1 
List of Frequencies 



ID 



Frequency 
(^Hz) 



Period 

(s) 



Amplitude 
(mma) 



Phase 
(s) 



fi> 

"15 

u 

ft 

"12 

/l 

"II 

h 
h 

"14 

A 
fa 
ft 
h 

"16 



1767. 
1774. 
1810. 
1894. 
1908. 
1933. 
1946. 
2059. 
2142. 
2214. 
2215. 
2328. 
2392. 
2521. 
2523. 
2659. 



3 ±0.3 
5 ±0.6 

7 ±0.2 

8 ±0.4 
1 ±0.7 

8 + 0.1 

9 ±0.5 
3 ±0.1 
8 ±0.3 
1 ± 0.3 
2±1.0 

8 ±0.2 

9 ±0.4 
7 ±0.4 
3 ± 0.3 
5 ±0.7 



565.84 ±0.09 
563.53 ±0.19 
552.26 ±0.06 
527.75 ±0.11 
524.08 ±0.18 
517.12 ±0.04 

513.65 ±0.14 
485.60 ±0.03 
466.67 ± 0.06 

451.66 ±0.05 
451.42 ±0.19 
429.41 ± 0.04 
417.91 ±0.07 
396.56 ±0.06 
396.30 ± 0.05 
376.02 ±0.10 



2.18±0.38 
1.28 ±0.43 
3.91 ±0.42 
3.07 ±0.87 
1.84 ±0.71 
8.51 ±0.82 
2.11 ±1.03 
4.59 ±0.36 
2.11 ±0.36 
5.44±1.15 
1.45 ±1.18 
2.88 ±0.35 
1.51 ±0.35 
3. 16 ±0.99 
3.86 ±0.99 
0.75 ±0.35 



194 ± 27 
12 ±56 
123 ± 18 
246 ±40 
344 ± 88 
139 ± 19 
466 ± 54 
82 ± 11 
97 ±22 
70 ±29 
31 ±107 
103 ± 15 
405 ±28 
31 ±37 
0±30 
26 ±50 



Notes. — Phases are times of maxima 
Frequencies with large uncertainties 



since HJD 2455495.5. 
are identified with u. 



squares fitting procedure, and all 16 frequencies, amplitudes 
and phases are left free in the final step. Several frequencies 
close to the most significant peaks do not converge to reliable 
values, and we consider these as uncertain (IDs labeled u in 
Table [TJ. It is not clear whether this is due to the poor cov- 
erage of our data, or caused by real amplitude variability of 
the target. Concern can also be raised about some frequencies 
that are barely resolved (e.g. fs and fa). 

5. DISCUSSION AND CONCLUSIONS 

A new, bright sdBV with an exceptionally low gravity 
has been found in a sample of UV-excess stars identified 
from GALEX photometry. In Figure |4] we show its position 
with respect to V338 Ser and the bulk of the sdB popula- 
tion, which falls between the zero and terminal age EHB 
lines. J20136+0928 is located close to minimum \ogg dur- 
ing the post-EHB stage of the uppermost evolutionary track, 
which corresponds to a 0.470 M helium burning star with 
an H-envelope mass of 0.004 M . Evolutionary models of 



more massive EH B stars are not shown, but can be found in 
lHan et alJ (120021) . from which one may infer that models with 
a mass of ~O.7M and an envelope of ~0.01 M* would ap- 
proximate the position of J20 136+0928 reasonably well. As 
mentioned in the introduction, such an overmassive EHB star 
cannot form through the regular binary mass-transfer scenar- 
ios that are responsible for producing the bulk of the EHB 
stars, but can form through the merger of two helium core 
WDs. Theoretical models for these scenarios, combined with 
future observational efforts, may allow us to determine the 
origins of both V338 Ser and J20136+0928. 

This discovery brings the total number of short period sdB V 
stars (both V361Hya a nd hybrid DWLyn st ars) up to 52 
(49 are summarised in I0stensen et al.1 1201 Obi one was re- 
cently discove red in the field of the K epler spacec raft by 
iKawaler et all 120101 and the last one by iBaran et "all 1201 lh . 
The discovery of pulsations in a star with such a low gravity 
as J20 136+0928 is well explai ned by the classical k mecha- 
nism for pulsations in sdBVs ( Charpinet et al.ll2001l) . How- 
ever, the exceptionally low gravity of J20 136+0928 is sur- 
prising in evolutionary terms, and the presence of pulsations 
allows the mass and internal structure to be explored. The en- 
couraging brightness of J20 136+0928 makes it an excellent 
target for time-resolved spectroscopy studies, and arranging a 
multi-site campaign with small telescopes when the target is 
visible throughout the night will be a priority for the coming 
season. By accurately determining the frequencies present in 
the star, asteroseismology will be able to reliably determine 
its mass and internal structure, and allow us to determine the 
most likely evolutionary origins of this object. 
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Based on observations made with the Mercator Telescope, 
operated on the island of La Palma by the Flemish Commu- 
nity, at the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astroffsica de Canarias. 
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